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(IT2FLC) is used. The controller consists of PD-Like FLC with integrated 
output. Kernek Mendel (KM) algorithm is used as the type reduction method 
Keyword: for the IT2FLC. This controller is applied to control a continuous binary 
ae trays distillation column. The controller has three tunable gains to reach 
Distillation nae ae me l 
minimum overshoot, minimum error and minimum settling time at least 
FLC variables can be controlled. The controller is a variable of the molar fraction 
IT2 of distillate and the reflex ratio is the manipulated variable. Integral Time 
Absolute Error (ITAE) is employed as an objective function to measure 
the improvement in time response where the error is between desired and 
output product composition. The performance of IT2FLC is compared with 
Type-1 PI Like FLC (TIFLC). The results of the simulations have shown 
that the project of IT2FLC works efficiently to no- disturbance and 
the effects of disturbance. Improve average is of 85% for a constant set-point 
without a disturbance and 80% with a disturbance. Furthermore, the average 
improvement for a step set-point is 53% without disturbance and 74% with 
disturbance. All results of the simulation confirmed the hardiness and control 
any consistent inaccurate with obvious advantages for the IT2FLC. 
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1. INTRODUCTION 

The columns of the distillation binary are likely the most of the controls in petroleum and chemical 
engineering. The dynamics of the variables thermo dynamic non ideality; size, uncertainty, and delay 
the measures of the composition put the control of procedures for the distillation with a challenging problems 
[1]. The auditors reported on fuzzy logic Type-1 that can be used successfully for a number of processes of 
non-linear and complex. Despite its wide popularity, it is shown in many scientific research cases the projects 
of the Type-1 FLC (TIFLC) have difficulties to minimize the impact of the uncertainties in the plant [2]. 

Therefore, the choice of the T1FLC does not necessarily the most favorable solution to the problem 
of control, whereas the best alternative is using Type-2 FLC (T2FLC) put the use of its fuzzy sets 
characterized by a consequence and, a larger number of parameters. T2FLC has been applied in many fields 
like vehicle nonlinear, suspension systems, robot mobile autonomous control, parallel mechanism of 
cable-driven, biochemical reactor control and inverted pendulum model. 

Pratishthananda and Worapradya [3] was developed the control with a real-time for a column of 
a binary distillation using Type-1 PI Fuzzy Logic Controllers. A nonlinear control adaptive algorithm for 
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a multi component batch distillation columns were developed by Jana and Murlidhar [4]. Canete et al. [5] 
applied a neurofuzzy adaptive control using soft sensors for continuous distillation column. Javadi and 
Hosseini [6] proposed the PI controller for a binary distillation column. Barceló-Rico et al. [7] realized 
the modeling and the control of a continuous distillation tower through fuzzy technique. A fuzzy model 
predictive control in multivariable control of distillation column was used by Sivakumar et al. [8] 
Al-Dunainawi and Abbod [9] performed three techniques of evolutionary algorithm as a tuner to fuzzy 
logic controller. 

Distillation column, as previously described, is a high nonlinear, poorly-defined relationship 
between the inputs and outputs, characterized by uncertainty. Thus; using T1FLC did not achieve high 
accuracy in controlling a binary trays distillation column and faced difficulties in reducing the effect of 
disturbances and uncertainties in the distillation column continuous type. 

The aim of this paper is to design IT2 PI-Like FLC for continuous distillation column (IT2FLC) to 
reach the set-point value for a distillate fraction (Xd) with minimum overshoot, minimum error and minimum 
settling time despite the disturbance. Results will be compared with using Type-1 PI-Like FLC taking into 
account the variations in system parameters and disturbance. 


2. DISTILLATION COLUMN PROCESS AND MODEL 

It is known that distillation columns are the units that most used in oil refineries; chemical and 
petrochemical plants [9]. It is the unit the operation in automatic control must be studied. There are a lot of 
types of distillation columns like, continuous or batch, binary or multi component, packed or tray [1]. 


2.1. Process description 

The distillation is a process in a liquid or vapor mixture of two or more components was separated 
into its fractions of required purity, by the application and removal of heat. Distillation column generally 
made up of trays of column. The distillation column consists of many trays, reboiler, condenser and a vertical 
column for the separation purpose. 

The feed is the entry in the column. Distillation column is separated into two sections. They are 
rectification section and stripping section. Trays under the feed tray are called the section of stripping and 
trays above feed tray are called rectification section [10], Figure 1 Leaving the Condenser liquid is entered 
into the drum of the reflex. A part of this liquid column deals as a reflex and fall down by the trays. Finally, 
the fluid is withdrawn, drum as a purified product. 
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Figure 1. Schematic impersonation of a typical binary continuous distillation column [11] 


2.2. Mathematical model 
The nonlinear system is modeled with the mass balance and component balance equations. 
The model is simplified under the following assumptions [12]: 
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— RRR /NG Chemical reactions occur. 

— No change in (constant) pressure. 

— Double (Binary) mixture. 

— No change in a (constant) relative volatility. 

— No vapor hold-up occurs in the trays. 

— Perfect equilibrium and mixing for vapor-liquid on all stages. 

— Perfect control level on reboiler and reflux drum. 

The Mass balance equation for each stage 1s [1]: 

— Consider stage 1=N + 1, xp=xn +1, where N is Number of stages, xp is the distillate liquid molar fraction, 
Xn is Liquid concentration on n‘” stage is a Condenser stage and is represented as (1): 


Mp——=V yi-1— L xp — D Xp (1) 


where Mp is Condenser hold-up, V is Vapor flow rate, yijis Vapor molar fraction on i" stage, L is Liquid flow 

rate and D is distillate flow rate. 

— Stage i=NF + 1 to i=N: represents trays above feed tray where NF is the number of the feed tray and its 
equation is as follows (2): 


Mi ZL xis V ya-p-L Xi- V yi (2) 


where M;is Hold-up on it” stage and X;is Liquid molar fraction on it” stage. 
— Stage i=NF: represents feed tray and its equation is written as follows (3): 


Mi “SL Xi +V ya-n-L XV yit F X; (3) 
where F is Feed flow rate and X¢is Feed liquid molar fraction. 
— Stages 1=2 to i=NF — 1: represent the stages below the feed tray which are expressed as follows (4): 


dxi 


Mi —=(L + F)xi+ı + V ya-p— L (L + F)xi— V yi (4) 


— Stage i=1, xp=x;: represents reboiler stage which is represented as follows (5): 


Ms “"=(L + P)xi+ ı + V yp- B xp (5) 


where Mg is Reboiler hold-up, xg is Bottoms liquid molar fraction, B is Bottom flow rate, and yg is Vapor 
molar fraction on bottom stage. 
— Vapor — liquid equilibrium relationship for each tray was represented as (6): 


a xi 

Yaa ©) 
where @ is Relative volatility 
— Reflux ratio is represented as (7): 

L 

where R is Reflex ratio. 
— Total reboiler mass balance is represented as (8): 

B=F- D (8) 
— Total condenser mass balance is represented as (9): 

V 
DSR a 
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2.3. Modeling of the distillation column 

The equations of the model are encrypted by the MATLAB environment, s-script function is called 
by the block in a Simulink. The resolution ordinary differential equations are used ODE45. Figure 2 shows 
the model Simulink of continuous type distillation column from (1 to 9). 
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Figure 2. Model of continuous type distillation column (1-9) 


Moreover, simulation model of the Distillation Column block is shown in Figure 3. 
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Figure 3. Simulink model of distillation column block (1-5) 


where the block CONSTILLss in the Simulink of distillation column is s-function that contains the program 
which describes equations (1 to 5). 


2.3.1. Simulation response of open loop distillation column 

The simulation was performed, to ensure compliance with the dynamic working in the operation, 
whereby the parameters are listed in the Table 1. The model has a non-linear reaction (response). Figure 4 
shows dynamic response with the change in the reflex ratio (RR) as +10% and -10%, where the nominal 
value of RR is 2.56 [1]. 
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Table 1. Parameters of distillation column model [1] 
Symbol Typical value and unit Symbol values and unit 


Mp 0.5 kmol F 1 kmol/h 
Mp 0.5 kmol a 2 

M; 0.5 kmol RR 2.56 
V 1.78 kmol/h N 20 

L 1.28 kmol/h NF 11 

D 0.5 kmol/h Xr 0.5 

B 0.5 kmol/h Xp 0.98 
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Figure 4. Dynamic response to the change in the reflex ratio (RR) 
(1) Response in the nominal value (Blue line), (41) Response at +10% in R(Red line), 
(111) Response at -10% in R (Yellow line) 


3. DESIGN OF THE PROPOSED CONTROLLER 
In the following section, the details of designing IT2 PI-Like FLC are explained. 


3.1. Design of IT2 PI-Like FLC (IT2FLC) 

The concept of T2FLC was proposed by Zadeh at 1975 as a statement of Type-1 FLC. Fuzzy sets 
function of the memberships are included in the Footprint Of Uncertainty (FOU), as shown in Figure 5. 
They are capable of modeling and handling the uncertainties, nonlinearities and linguistic associated with 
the inputs and output of FLC by modeling them and minimizing their effect. Therefore, IT2FLC is used 
instead of the TIFLC. IT2 fuzzy sets are extensions of the classical fuzzy sets that add an extra degree of 
freedom to entitle an uncertainty factor in the system but they are not easy compared with Type-1 FLC [13]. 

In this paper, IT2 PI like FLC of Mamdani type is designed with seven inputs-output triangular 
shaped memberships. The universe of discourse range is within [-1 to 1] for each input and output. 
The forms of MFs are shown in Figure 6. Centroid is the method that is selected for the defuzzification 
technique. The Rule-base of a proposed IT2FLC is shown in Table 2. KM algorithm is used to implement 
the type reduction method of the IT2FLC. The inputs and output gains of the controller with the rules were 
tuned to reach the best response to reach minimum overshoot, minimum settling time and minimum error in 
distillate molar fraction. The features of using Mamdani type are: it is my thought, it suits human input and it 
is commonly accepted. 
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Figure 5. IT2 membership function structure [14] Figure 6. Input/Output MFs 
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Table 2. Rule base of IT2FLC 
e/Ae NB NM NS Z PS PM PB 


A PI-like IT2FLC probably employed to serve as a PD such as IT2FLC. As shown in (10) shows 
a PD controller equation results in a position form, however (11) shows a PI controller [15]. 


Up(k) = KPe(k)+KDAe(k) (10) 
RR(k) = Kout.(ulo(k)+uD(k)) (11) 


where e(k) is a signal of sampled error, Ae(k) is the sampled error signal rate change, index k represents 
the present sampling instant, RR(k) is a control action, Kout is an output gain, uJo(k) is the previous an 
control action and up(k) is derivative control action. 


3.2. Design of the controller logic inaccurate type 1 (TIFLC) 

This article describes a Mamdani T1 PI Like FLC [1] is proposed as a comparison with IT2FLC 
where it should with two inputs and one output. The entries are going through a process of fuzzification with 
(5) membership. The trapezoidal membership was selected in the extremes and triangular membership in 
the middle and the defuzzification selected as the method is Centroid. More details can be found in 
the reference 1. 


4. RESULTS AND DISCUSSION 
The controlled system is tested according to the following scenarios: 
— Constant Distillate fraction set-point (0.98) 
— Constant Distillate fraction set-point (0.98) with disturbance 
— Step Distillate fraction Set-point 
— Step Distillate fraction Set-point with disturbance 
The Distillate fraction (XD) set-point is selected from a Table 1 as a purity of 98%, the time of 
the simulation is 250 hours and the function of the cost used for the evaluation of the performance of 
the control is ITAE: 


J=250 f 


le|dt (12) 

The applied disturbance is the Feed flow rate (F) and Feed liquid molar fraction XF [1], as shown in 
the Figure 7. The feed flow rate is represented as a Sine wave oscillation with low frequency and the Feed 
liquid molar fraction is represented as a smooth square wave with high frequency. 





Tirri 





(b) Feed liquid molar fraction 


Figure 7. Disturbances on the entire feed of the distillation column: 
(a) Feed flow rate (b) Feed liquid molar fraction 
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The simulation of a closed loop for a distillation process using IT2FLC is shown in Figure 8 and 
the simulation of the closed loop for a distillation column process using IT1FLC is shown in Figure 9. 






























































Figure 9. Closed loop system for a distillation column process using IT1FLC 


The gains of the proportional control (KP), Derivative control (KD) and the output gain (Kout) are 
tuned manually for both controllers to reach the best time response as shown in Table 3 then the simulation 
response for the Xd concentration in both constant set-point and step set-point with disturbance and without 
disturbance for a Type-1 FLC are shown in Figure 10 and Figure 11 respectively and for an IT2FLC are 
shown in Figure 12 and Figure 13 respectively. 
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(a) Step set-point only (b) Step set-point with disturbance 


Figure 11. Time trajectory response of the system Xd consideration with fuzzy Type-1 for 
(a) Step Set-point only (b) Step Set-point with disturbance 
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Table 3. Tuned gains for IT2LC 


Types PI-Like IT2FLC PI-Like IT2FLC PI-Like IT2FLC PI-Like IT2FLC 
With constant set- With constant set-point and With step set- With step set-point 
Gains point only disturbance point only and disturbance 
KP 5 17.5 5 8.179 
KD 5 9.98 5 3.4197 
Ku 10 9.8 30 8.9 
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Figure 12. Time trajectory response of the system Xd consideration with fuzzy Type-2 for (a) Constant set- 
point only (b) Constant set-point with disturbance 
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Figure 13. Time trajectory response of the system Xd consideration with fuzzy Type-2 
(a) Step set-point only (b) Step set-point with disturbance 


By applying 12, the percentage of enhancement between the IT2FLC compared with ITIFLC for 
different cases are listed in Table 4. 


Table 4. Percentage of enhancement between the IT2FLC compared with TIFLC 


Contoiler Constant Step Set- percentage of enhan-cement Percentage of 
Set-point point for constant Set-point enhancement Step input 
Type-2 76.23 286 ~85% ~53% 
Type-2 with disturbance 141.1 470.6 ~80% ~7T4% 
Type-1 506.8 604 
Type-1 with disturbance 131.3 1831 
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5. CONCLUSION 

The application of IT2 PI-Like FLC to the nonlinear dynamical model of a continuous type 
distillation column demonstrates the effectiveness in reaching a top control performance and allowing a faster 
and exact control of the process. The disturbance rejection and set point tracking is done with less amount of 
overshoot as a compared with using Type-1 PI-Like FLC. Simulations carried out, showed an average 
improvement measured using ITAE of 85% for a constant set-point without disturbance and 80% with 
a disturbance, for a step set-point the improvement 1s 53% without disturbance and 74% with disturbance. 
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